Exactly quantized dynamics of classical incommensurate sliders 
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We report peculiar velocity quantization phenomena in the classical motion of an idealized ID 
solid lubricant, consisting of a harmonic chain interposed between two periodic sliders. The ratio 
Vcm/vext of the chain center-of-mass velocity to the externally imposed relative velocity of the sliders 
stays pinned to exact "plateau" values for wide ranges of parameters, such as sliders corrugation 
amplitudes, external velocity, chain stiffness and dissipation, and is strictly determined by the 
commensurability ratios alone. The phenomenon is explained by one slider rigidly dragging the 
kinks that the chain forms with the other slider. Possible consequences of these results for some real 
systems are discussed. 

PACS numbers: 68.35.Af, 05.45.Yv, 62.25.+g, 62.20.Qp, 81.40.Pq, 46.55.+d 



In this Letter we report surprising exact velocity quan- 
tization phenomena in a one-dimensional (ID) non-linear 
model inspired by the tribological problem of two slid- 
ing surfaces with a thin solid lubricant layer in between 
[ij, |6j. The model layer consists of a chain of N har- 
monically interacting particles interposed between two 
rigid generally (but not necessarily) incommensurate si- 
nusoidal substrates (the two "sliding crystals" , sketched 
in Fig. ^) externally driven at a constant relative ve- 
locity Vf,xt. The equation of motion of the i-th particle 



mxi 



sinfc+(xi — v+t) -f F_ sink_{xi — v_t)\ 
± 



2 

-I- K{x.i+i -I- Xj-i - 2x^) - 7 ^ - v±) 



(1) 



where m is the mass of the N particles, K is the chain 
spring constant, and k± — (27r)/a± are the wave- vector 
periodicities of potentials representing the two sliders, 
moving at velocities v±. We set, in full generality, v+ = 
and Vf-xt = u- — u+ . 7 is a phenomenological parameter 
substituting for various sources of dissipation, required 
to achieve a stationary state, but otherwise with no ma- 
jor role in the following. F± are the amplitudes of the 
forces due to the sinusoidal corrugation of the two sliders 
(we will commonly use F_/F^ = 1 but we checked that 
our results are more general). We take a+ — 1, m — 1, 
and = 1 as our basic units. The relevant length-ratios 
0,0] are therefore r± = a±/ao; we will take, without loss 
of generality, r_ > r+ , and confine our attention mostly 
to cases with r+ > 1. The inter-particle equilibrium 



length oo, not entering explicitly the equation of motion 
appears only via the boundary conditions, which are 
taken to be periodic (PBC) 0, XN-i-i = xi -\- N qq, to 
enforce a fixed density condition for the chain , with a 
coverage r-)_ of chain atoms on the denser substrate. Pre- 
vious studies on a related model |3| achieved the sliding 
through the application of a constant driving to one of 
the two substrates, via an additional spring. That pro- 
cedure obscures the surprising quantization phenomena 
which are instead uncovered when sliding occurs with a 
constant velocity Vg^t- 

Upon sliding the substrates, Vext 7^ 0, the lubricant 
chain slides too. However, it generally does so in a very 
unexpected manner: the time-averaged chain velocity 
w — Vcm/vext, IS generally asymmetric, namely different 
from 1/2. Even more surprisingly, w is exactly quantized, 
for large parameter intervals, to plateau values that de- 
pend solely on the chosen commensurability ratios. The 
asymmetrical w-plateaus are generally very stable, and 
insensitive to many details of the model, due to their in- 
trinsically topological nature; we show that they are the 
manifestation of a certain density of solitons in the lubri- 
cant which are set into motion by the external driving. 

We now turn to illustrating these results. We inte- 
grated the equations of motion (Q) starting from fully 
relaxed springs (x^ = iag, ij = v^xt/^), by a standard 
fourth-order Runge-Kutta method. After an initial tran- 
sient, the system reaches its dynamical stationary state, 
at least so long as 7 is not exactly zero. Figure^ shows 
the resulting time-averaged center-of-mass (CM) velocity 
Vcm as a function of the chain stiffness K for four repre- 
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FIG. 1: (Color online) (a) Average drift velocity ratio 
w — Vam/vext of the chain as a function of its spring stiff- 
ness K for different length ratios (r+,r_): commensurate 
(3/2,9/4), golden mean (GM) (</>, 9!)^), spiral mean (SM) 
(a, (T^) (a ~ 1.3247. . . is the solution of cr'^ = a + 1), and 
The plateau labeling is explained in the text. Here 
7 = 0.1, Vext ~ 0.1, and periodic boundary conditions 01 are 
used. The (0, 0^) 1/1 plateau value is w — 0.381966 . . . , iden- 
tical to 1 — to eight decimal places, (b) The main plateau 
speed w as a function of r+. (c) A sketch of the model. 

sentative (r+,r_) values. We find that w is generally a 
complicated function of K, with flat plateaus and regimes 
of continuous evolution, in a way which is qualitatively 
similar for different cases. The main surprise is that all 
plateaus show perfectly flat w = Vcm/vext that are con- 
stant (quantized) to all figures of numerical accuracy, the 
precise value strikingly independent not only of K, but 
also of 7, Vext , and even of i<L / . Open-boundary simu- 
lations show moreover that the PBC used are not crucial 
to the plateau quantization, which occurs even for a lu- 
bricant of finite size and not particularly large iV, such 
as a hydrocarbon chain molecule would be. 

To explore the origin of the universality of w, we ana- 
lyze the dynamics for a large number of values of (r+ ,r^), 
and observe that: (i) at least one velocity plateau as a 
function of K occurs for any (r_|_ , r_ ) ; (ii) additional nar- 
rower secondary plateaus often arise for stiffer lubricant 
(larger K, see Fig.QJ; (iii) the velocity ratio w of the first 
plateau found for increasing K, satisfies w = 1 — for 
a large range of (r+, r_). 
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FIG. 2: (Color online) Snapshots of the distance between 
neighbor lubricant particles in the chain Xi — Xi-i at three 
successive time frames. All parameters as in Fig. Q except 
r_ ~ 10.36, and A' = 10 (inside the main plateau), (a) 
r+ = 1.031 (kink density S/a+ = 0.031); (b) r+ = 0.995 
(anti-kink density |5|/a-i- = 0.005). 

We can understand these results as follows. Con- 
sider initially the situation of quasi-commensuration of 
the chain to the a+ substrate: — 1 + S, with small 
S. This induces a density psoi — S/a+ of solitons (or 
kinks, essentially substrate minima holding two parti- 
cles, rather than one) [sj. The second, less oscillat- 
ing a„ slider, which moves at velocity Ve^tj will drag 
the kinks along: Vgoi = Vext- If Po = 1/ao = r+/a+ 
is the linear density of lubricant particles, mass trans- 
port will obey VcmPo = VsoiPsoi- This yields precisely 
w = Vcm/vext = Psoi/po = S/r+ = 1 - Thus the 

exact plateaus arise because the smoother slider (whose 
exact periodicity r_ > r+ is irrelevant) drags the kinks, 
of given density, at its own full speed Vext , as illustrated 
Fig.Et- As shown in Fig.^, this physics extends even to 
large \S\ 1, where no individual kink can be singled out. 
This works even for 6 < (anti- kinks), where, remark- 
ably, the lubricant CM moves in the opposite direction 
with respect to the driving v^xt (w < 0, see Fig. ^ for 
r-|_ — 4)~^): exactly as holes in a semiconductor, anti- 
kinks (carrying a negative "charge") moving at velocity 
+Vext effectively produce a backward net lubricant mo- 
tion. The motion of the anti-kinks (regions of increased 
inter-particle separation) is illustrated in Fig. ^p. 

The motion of particles is instructive. Figure O plots 
the time evolution of the velocity of a single chain particle 
Xi, and of Vcm, for a value of K inside a plateau, K = 1, 
for r+ = (/) and r+ cr of Fig.^ Here two clear kinds of 
behavior emerge. Remarkably, the single-particle motion 
in the GM plateau is perfectly time-periodic. A simi- 
lar periodic dynamics is found, in appropriate regimes, 
for all rational and quadratic irrational (r+ , r_ ) values 
we tested. The SM motion is definitely not periodic. 
The Fourier spectrum of the particle motion, shown in 
Fig. 121;, confirm that the SM yields only quasi-periodic 
orbits with two prominent incommensurate frequencies 
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(c) Spiral Mean: v, FFT 
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FIG. 3: (Color online) Time evolution of a particle velocity 
Xi, and of the chain CM velocity Vcm (fluctuations rescaled by 
a factor 50), for the GM (a) and SM (b) cases of Fig.[Tl Ampli- 
tudes of the Fourier spectrum of Xi{t) (c) and of Vcm{t) (d) for 
r = a. Individual particle spectra have identical amplitudes, 
and differ only in the phases, which leads to a remarkable can- 
cellation in the Vcm power spectrum. Here K = 1, y = 0.1, 
and Vext = 0.1. 
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FIG. 4: (Color online) Dynamical depinning force Fc for 
the GM plateau (a), extracted by a slow adiabatic increase 
of the force F applied to all particles, until intermittencies 
appear, signaling collective slips (b). The critical behavior 
of the average velocity for K — > Ka (c,d), showing the typi- 
cal square-root singularity (dotted line) associated to type-I 
intermittencies. Here 7 — 0.1, and Vext = 0.1 were used. 



/+ and /_. However, a phase-cancellation between the 
Fourier spectra of different chain particles (aU having the 
same amplitude spectrum, with different phases) yields 
strictly periodic CM motion. Periodic CM oscillations 
around an exactly quantized drift velocity is a common 
feature of all plateaus in the chain dynamics. These peri- 
odic oscillations can be understood as the solitons moving 
at velocity v^xt encountering a periodic Peierls-Nabarro 
potential of period . The corresponding frequency 
of encounter Vext/a+ is clearly visible in Fig.|31i. We ob- 
serve that the Peierls-Nabarro barrier vanishes (strictly 
constant Vcm) for all cases where the particle motion is 
periodic. 

To understand why single particles may move peri- 
odically in some plateaus, let /+ = Wcm/a+ and /_ = 
{vext — Vcm)/oL- be the average encounter frequency of 
a generic particle with the two substrate periodic corru- 
gations. Whenever these two frequencies are mutually 
commensurate, 

n+/+=n_/_ (integers n±) , (2) 

each particle is driven by the two generally incommen- 
surate corrugations in an effectively periodic way, and 
undergoes a periodic oscillation, of period T = n^/ f+ = 
n+/ f-. This is observed in the Fig. ^ examples: r+ — cj) 
yields n^/ri- = 1/1 and = 3/2 yields n^/n^ = 4/3. 



The plateaus in Fig. are labeled accordingly. When 
/-!_ and /_ are incommensurate, the motion is quasi- 
periodic, as illustrated by the SM case of Fig. |31 By 
definition, /+//- = {r^ /r^) w/ {I — w): for any given 
plateau velocity w, suitable choices of r^/r^ can make 
the two frequencies either commensurate, as in Eq. jSJ, 
or incommensurate. In particular, for the main plateau 
we have /+//- = r_(r+ — l)/r-)_, which can be made 
rational for any choice of r+, by choosing a suitable 
r_ = r-|_(r+ — 1)^^ rt_/rt+. For example, we verified that, 
for r+ = cr and r_ = o'(o'— 1)^^, individual particles do 
oscillate periodically in the main plateau. 

Low driving velocities Vext are beneficial to the ap- 
pearance and width of plateaus. For increasing Vext, the 
plateaus shrink and eventually disappear, still remain- 
ing exact while they do so. The critical Vext where the 
plateaus end depends on K, but is usually smaller than 
unity, for the parameters of Fig. The 1/1 plateau of 
the GM case is especially wide and robust against an in- 
crease of Vext (for the parameters of Fig. and K — A, 
up to Vext — 1-5) and other perturbations. This seems 
related to (j) being uniquely associated to the equal drive- 
frequency ratio, n+/n^ = 1/1: cf) appears therefore, in 
the present dynamical context, as the "most commen- 
surate^'' irrational, at variance with static pinning in the 
standard FK model, where the opposite is true IJjQ- 

The finding of exact plateaus implies a kind of "dynam- 
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ical incompressibility" , namely identically null response 
to perturbations or fluctuations trying to deflect the CM 
velocity away from its quantized value. In order to probe 
the robustness of the plateau attractors, we introduced 
an additional constant force F acting on all particles in 
the chain, trying to alter the force-free sliding in the 
plateau regime. As expected, as long as F remains suffi- 
ciently small, it does perturb the single-particle motions 
but has no effect whatsoever on w, which remains exactly 
pinned to the attractor {vcm = Vpiateau)- The plateau dy- 
namics is only abandoned above a critical force Fc- This 
dynamical depinning takes place through a series of type- 
I intermittencies 0, as shown in Fig.^js where Vcm{t) is 
plotted against a slow adiabatic ramping of F. 

A precise value of Fc can be obtained by ramping F 
with time with a gentle enough rate of increase or, alter- 
natively, by a Floquet-Lyapunov linear stability analysis 
0, to be shown elsewhere. The value of Fc is a func- 
tion of the parameters, and Fc vanishes linearly when K 
approaches the border Kc of the plateau, as in Fig. 
The depinning transition line Fc, ending a,t K — Kc, is 
a 'first-order' line, with a jump Av in the average Vcm 
and a clear hysteretic behavior (not shown) as F crosses 
Fc- As can be expected, Av decreases to as X in- 
creases towards Kc- Thus K = Kc represents a genuine 
non-equilibrium critical point, where the sliding chain 
enters or leaves a dynamical orbit. The precise value of 
Kc depends on parameters such as Vext and 7; however, 
its properties do not. As K approaches Kc from above 
(no external force), Vcm approaches Vpiateau m a criti- 
cal manner, as suggested in Fig. This is detailed in 
Fig. EIc,d), where the critical behavior is shown to be 
Av oc {K — KcY^"^ , the value typical of intermittencies 
of type I 0. For K > Kc, in fact the chain spends most 
of its time moving at Vcm (t) ~ Vpiateau , except for short 
bursts at regular time-intervals t, where the system as a 
whole jumps ahead by Oq, i.e. an extra chain lattice spac- 
ing (collective slip). The characteristic time t between 
successive collective slips diverges as r cx (X — Kc)^^^^ 
for K Kc, consistent with the critical behavior of w. 
We verified that the if-plateaus for more general values of 
r+ and of r_ show the same kind of infinite stiffness, and 
a critical decrease of Fc near the plateau edge, similar to 
that of Fig. Hi for the GM. 

The phenomena just described for a model ID system 
are quite extraordinary; it would be fascinating if they 
could be observed in real systems. Nested carbon nan- 
otubes [inj |. or confined one-dimensional nanomechan- 
ical systems possible arena for the phe- 

nomena described in this Letter. Though speculative 
at this stage, one obvious question is what aspects of 
the phenomenology just described might survive in two- 
dimensions (2D), where tribological realizations, such 
as the sliding of two hard crystalline faces with, e.g., 
an interposed graphite flake, are conceivable. Our re- 
sults suggests that the lattice of discommensurations - a 



Moire pattern- formed by the flake on a substrate, could 
be dragged by the other sliding crystal face, in such a 
manner that the speed of the flake as a whole would 
be smaller, and quantized. This would amount to the 
slider "ironing" the solitons onward. Dienwiebel et al. 
|l2l | demonstrated how incommensurability may lead to 
virtually friction-free sliding in such a case, but no mea- 
sure was obtained for the flake relative sliding velocity. 
Real substrates are, unlike our model, not rigid, sub- 
ject to thermal expansion, etc. Nevertheless the ubiq- 
uity of plateaus shown in Fig. ^ and their topological 
origin, suggests that these effects would not remove the 
phenomenon. A real-life situation with a distribution 
of differently oriented crystalline micro-grains, each pos- 
sessing a different incommensurability, is also potentially 
interesting; each grain, we expect, will tend to stabilize a 
certain average CM velocity depending on its incommen- 
surability. Other realizations or applications inspired by 
the physics described by our model might be accessible, 
notably in grain boundary motion, in the sliding of opti- 
cal lattices 13j or of charge-density- wave systems 
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